Addition Spectra of Chaotic Quantum Dots: 
Interplay between Interactions and Geometry 
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We investigate the influence of interactions and geometry on ground states of clean chaotic quantum 
dots using the self-consistent Hartree-Fock method. We find two distinct regimes of interaction 
strength: While capacitive energy fluctuations 5x follow approximately a random matrix prediction 
for weak interactions, there is a crossover to a regime where Sx is strongly enhanced and scales 
roughly with interaction strength. This enhancement is related to the rearrangement of charges into 
ordered states near the quantum dot edge. This effect is non-universal depending on dot shape and 
size. It may provide additional insight into recent experiments on statistics of Coulomb blockade 
peak spacings. 
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The understanding of the interplay between many- 
body interactions and chaos has evolved to a promi- 
nent field in mesoscopic physics. This challenging is- 
sue brings together two seemingly disconnected fields, 
namely many-particle physics and quantum chaos. Clas- 
sically chaotic independent-particle dynamics is found in 
mesoscopics for both, disordered systems due to scat- 
tering at impurities and clean quantum dots, where the 
shape of the confining potential can give rise to chaotic 
electron motion. While interaction effects in disordered 
systems have been intensively studied in the recent past 
PI , accounts on the inter-relation between many-body ef- 
fects and chaotic dynamics in ballistic quantum dots are 
still rare (see, e.g., Refs. 1^,^). The latter systems can be 
modelled by quantum billiards which have served as pro- 
totypes to investigate quantum signatures of integrable 
and chaotic single-particle dynamics m. Hence a gener- 
alization of these quantum chaos concepts to interacting 
particles in billiards appears natural. In this letter we 
approach this problem within a self-consistent Hartree- 
Fock (SCHF) approximation after checking its validity 
by comparing with exact diagonalization. 

Our work was partly motivated by new measurements 
of charge transport through quantum dots with variable 
size and shape, which represent ideal tools to investigate 
experimentally interaction effects in confined chaotic sys- 
tems. If such a dot with TV electrons and ground state 
energy En is only weakly coupled to external leads a con- 
ductance peak is observed, upon tuning a gate voltage 
connected to the dot, whenever the chemical potential 
/i^r — En — En-\ coincides with that of the leads. In 
between the pronounced conductance peaks the current 
through the dot is suppressed due to Coulomb blockade 
[p|. The spacing between neighboring peaks as function 
of gate voltage is proportional to the capacitive energy 
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A number of recent experiments |q-pf| showed that fiuctu- 
ations of xn resemble a Gaussian distribution, while the 



assumption of a constant interaction (CI) , which properly 
describes the mean peak spacing, combined with random 
matrix theory (RMT), predicts a Wigner-Dyson distribu- 
tion. Two reasons have been proposed for the Gaussian 
distribution, the scrambling of the single-particle spec- 
trum due to interactions |6|Jl(]|-[l^ and the deformation 
of the dot by changing the gate voltage [Q . 

However, it is not yet understood why the measured 
peak spacing fluctuations, Sx = Hxn) ~ (xn)'^)^^'^, sig- 
nificantly vary between the different experiments: Pa- 
tel et al. ||] found rather small fluctuations which are 
comparable to the single-particle mean level spacing A. 
They are in line with the RMT-fCI model predicting 
Sx = A(4/7r - 1)1/2 ^ 0.52A, and with RPA calcu- 
lations for weak interactions (high densities, r^ < 1) 
|la,E6[. On the contrary, the other experiments yield 
large fluctuations which scale rather with (xn) than with 
A: Sx ~ 0.06 - 0.15(xjv) in Refs. §0§. E.g., in the 
recent experiment B , where the shape of the dot is ex- 
pected not to be distorted by tuning the gate voltage, Sx 
was found to be 15 times larger than predicted by RMT. 

Thus two main open questions arise: (i) what are the 
responsible mechanisms for the large fluctuations and (ii) 
why do they show such a sample-dependent behavior, 
assuming that the samples are chaotic? 

Here we consider clean chaotic quantum dots and study 
the interplay between the Coulomb interaction of the 
electrons and the geometry of the dots. We propose a 
mechanism for enhanced capacitive energy fluctuations 
Sx in chaotic structures and show that: 

(i) interaction drives the systems to a regime of large, 
non-universal fluctuations by developing ordered charge 
states near the edges. 

(ii) the shape of the confinement is an important fac- 
tor for the formation of such states and hence for the 
fluctuations. 

To our knowledge, all related theoretical studies con- 
sider interaction effects for disordered models P,nG-n3|, 



apart from the work by Stopa |g], who suggested the 
occurence of strongly scarred states as a mechanism for 
enhanced fluctuations. Our calculations show that 5x 
can be enhanced in the absence of scarred states. 

We model the quantum dots by a family of billiards 
which arise from a deformation of the disk |]l^. They 
have already been employed as non-interacting models 
in the context of Coulomb blockade |14 Iq]. Here they 



will serve as appropriate tools to study interaction ef- 
fects upon changing the billiard geometry. The con- 
finement is defined by a hardwall potential U satisfy- 
ing U{u,v) = inside a domain D and U{u,v) — oo 
outside. Using complex coordinates lo — u + iv and 
z = X + iy, D is defined by the conformal mapping |17| ] 
uj{z) = R{z + bz^ + ce^^z^)/{Vl + 2b^ + 3c'^); |z| < 1 . 
The real parameters 6, c, and S determine the shape, and 
R defines the dot size and area, A = ttR^ . For the de- 
formations to be considered (insets in Fig. ^) the classi- 
cal single-particle dynamics is predominantly chaotic (de- 
scribed by a large Kolmogorov entropy for the strongest 
deformation [Q), and we found that the single-particle 
energies exhibit Wigner-Dyson statistics. 

To include interactions we first calculate the single- 
particle eigenfunctions <j)i and -energies e^ of a de- 
formed billiard. We then use them to construct 
the Coulomb interaction matrix elements Viju = 
j dridr2(t>i{vi)4>j{v2)V{vi -V2)(t>k{ri)(j)i{r2) with V{t) = 
e^/er, and the many-electron Hamiltonian 
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Here cl{ci) creates (annihilates) the i-th eigenstate of the 
non-interacting Hamiltonian. In terms of energy units 
h?/{2m*R'^), Vijki is proportional to the dimensionless 
interaction strength, or system size, R/a^, where a^ = 
h^€/m*e'^ is the effective Bohr radius. Eq. (||) can serve 
as a starting point for both exact diagonalization for few 
electrons |1£| and SCHF calculations. 

In this work we focus on ground state properties for 
many-electron quantum dots and employ SCHF system- 
atically for systems up to 27 spinless electrons. The 
SCHF allows us to calculate the ground state energies 
even in the presence of two-electron density correlations 
which are ignored in the RPA perturbation calculations 
0,[6|. To check the validity of our SCHF results for 
the capacitive energies xn we first compare them with 
exact diagonalization results for few electrons. As shown 
in Fig. n the SCHF results are very accurate even for 
strong interactions. 

The SCHF results for xn do not show any visible de- 
pendence on N but rather on the interaction strength 
(system size) R/a*g. Let us thus consider the mean fluc- 
tuations Sx, where the statistical average (. . .) is per- 
formed over N as in the experiments 



Our main results are summarized in Fig. g. It shows 
dx as a function of interaction strength for three dif- 
ferent geometries depicted as insets close to the curves. 
For small size quantum dots, R/a*^ ^ 12, 5x is close 
to the RMT prediction. Most interestingly, we find for 
larger R/a*g a cross-over to a regime where the capac- 
itive energy fluctuations increase roughly linearly with 
increasing interaction strength. To our knowledge such a 
clear-cut cross-over behavior has not been found in any 
tight-binding model for disordered systems. This indi- 
cates that it could be a property of the charge distribu- 
tion in confined ballistic systems. 

Using the standard definition of r^, namely r^ = 



A/{TTNa*g'^), we have rg = {R/a*g)/\/N. Since we cal- 
culate 5x over a range 3 < A'' < 27 for a given value of 
R/a*g in Fig. 0^s has not a well-defined value but is in 
the range {l/V27)R/a*g < r^ < (l/V3)i?/a|j. Hence the 
cross-over value of R/a*^ in Fig. implies that 5x retains 
its RMT value up to Tg « 2 which is consistent with the 
RPA approach of Ref. ||l|l. 

In the following we discuss the mechanism responsible 
for the cross-over to enhanced fluctuations 6x- To this 
end it proves convenient to decompose xn+i as pq] 



Xn+1' 



XN+l) (N) 



'N+2 



^N+2 



-JV+2 ^iV+l 



= Ei+E2. (3) 
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Here, ej denotes the j-th eigenvalue of the SCHF Hamil- 
tonian with i electrons. Quantitatively, Eq. y^, which 
relies on Koopmans' approximation, does not give the 
same value as 8x in Fig. g. However, Eq. (||) shows the 
same qualitative behavior with the same cross-over point 
but slightly larger slope than in the SCHF case. Fig. g. 




FIG. 1. Comparison of self-consistent Hartree-Fock 

(SCHF) and exact diagonalization (EXACT) results of ca- 
pacitive energies xn, Eq. (|l|), as a function of interaction 
strength for A'^ = 3 (upper) and N = A (lower curve) for the 
billiard geometry shown in the bottom inset of Fig. S. 

The quantity Ei in Eq. (0), which describes the in- 
teraction energy between the electrons in (A'^-t-l)th and 
(A + 2)th state, has been considered in Ref. ||l^. £'2 is 
the energy difference between two adjacent single-particle 
levels of the same Hamiltonian. In the RPA approach 
|l5| , it was claimed that {E2) is given by A and its fluctu- 
ations 5E2 follow the RMT value ~ 0.52A. As displayed 
in the inset of Fig. ^ our calculations show that this seems 
to be true for R/a*^ ^ 12, which includes the regime 



Vs < 1, where RPA is reliable. However, {E2) / ^ grows as 
the system size increases, and 6E2/A and 5-^2/(^2) do 
no longer follow the RMT prediction for R/a*g > 12. 

The similarity of the cross-over behavior of 5x and dE2 
in Fig. ^ suggests that an analysis of E2 proves appropri- 
ate to understand the enhanced fluctuations Sx- 

What happens in the regime R/a*g > 12? To answer 
this question we present in Fig. the evolution of unoccu- 
pied SCHF single-particle levels ei^\n = N+l, N+2, . . . 
with increasing interaction for a few representative val- 
ues of N. We find that the enhancement of mean and 
rms values of E2 = {^^^2 ~ ^jv-i-i) i^ accompanied by the 
appearance of sudden jumps of the single particle levels 
as marked by vertical arrows in Fig. ^ (Note however 
that not for all values of N the energy levels show such 
jumps in the R/a*^ range considered, e.g. A^=16.) 
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FIG. 2. Fluctuations of capacitive energies Sx ~ Hxn) 
N)^y as a function of interaction strength for inter- 
acting electrons in three chaotic billiards of different shape 
shown on the right (defined with 5 = 1.5, fe = c = 0.12, 0.16, 
and 0.2 from above, see the text). 5x shows a clear cross-over 
near R/a*g — 12 and increases with increasing R/ag in the 
strong interaction regime. An analysis of the statistical error 
in computing 5x shows that it is of the order of the fluctu- 
ations visible in the curves. Inset: The mean and rms value 



of £2 = e 
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for the quantum dot of the bottom right 



inset. Note that for larger interactions {E2) 7^ A. 

The abrupt changes in the SCHF single-particle ener- 
gies go along with a rearrangement of the corresponding 
wavefunctions as shown in Fig. 0(a,b). This indicates a 
change in the total charge densities and thereby in the un- 
derlying mean-field potential. The two lower left (right) 
panels in Fig. depict the total charge densities for in- 
teraction strength R/a*g = 5 (15) for two selected, rep- 
resentative N. A systematic analysis of a large number 
of geometries and particle numbers allowed us to extract 
two common trends visualized in the figure: 

(i) With growing system size the electrons tend to be 
localized near the edge of the dot. For R/a*g > 12 the 
charges are reorganized in an ordered structure similar to 
a one- dimensional crystal or a charge density wave state. 



(ii) An abrupt change in the mean-field potential for 
fixed N occurs when an electron from the inner region of 
the dot moves into a state close to the edge. 

Let us discuss the connection to the increase in {E2) 
and 6E2- Electrons in states extending over the whole 
dot for small R/clq move with increasing interaction 
to the edge due to the long-range Coulomb repulsion. 
When the states near the edge are fully occupied (see 
e.g. Fig. ^(d)) the corresponding mean field reduces the 
accessible space for the next electron above the Fermi 
level (crossover from ^(a) to (b)). The reduced effec- 
tive size of the dot gives rise to an enhanced level spac- 
ing {E2) which is accompanied by enhanced fluctuations. 
Note that 6E2/{E2) also deviates from the RMT value of 
« 0.52. This indicates that even for a chaotic geometry 
the SCHF single-particle states near the Fermi level in 
the strong interaction regime do no longer follow RMT 
due to charge rearrangement near the edge. 

The reorganized charge states at the edge differ from a 
two-dimensional Wigner crystal where the electrons are 
localized due to electron correlation energy Ig^l ■ 

One can distinguish two kinds of ordered states: 

For .smaller N a phase which might be called a "crystal 
chain" evolves where the electrons are localized near the 
edge with clear spatial separation (see, e.g.. Fig. 0(d)). 
Due to quantum fluctuations of such a ID crystal lattice 
this phase will be bounded from above by some critical 
electron number Nc- Since the minimum distance be- 
tween the electrons should be larger than a Fermi wave- 
length Xp , the maximum number of electrons which can 
reside at the edge is « Lp/Xp, where Lp is the perimeter 
of the dot. Since in the crystal chain state all electrons 
are at the edge, Nc can be estimated as Nc ~ Lp/A, 
where we used Xp ~ R/VN. For the dot geometries 
used we find Nc to be of order of 10. This rough estima- 
tion is in line with our numerical results where for N — 7 
all electrons are near the edge, whereas for A^ = 11 only 
10 electrons seem close to the edge with one electron near 
the center of the dot (see Fig. |j(f)). 

For /ar(7er electron number N>Nc, as in Fig. Q(f), the 
charge density exhibits a modulation close to the edge 
rather than well separated charges. Such ordered, charge 
density wave like structures [£1| persist up to high elec- 
tron numbers. We still found these states for A^ = 40, the 
highest N considered. They play the same role for the 
enhancement of 5x (see, e.g.. Fig. y for A^ = 24) as the 
crystal type states do for N<Nc. 

A few further remarks are due: 

(i) A strongly deformed quantum dot where part of the 
boundary is concave [E2| (insets in Fig. 0) renders the 
formation of edge states more difficult than the dot with 
a shape close to a disk due to the interplay between the 
long range interaction and geometry. This is the reason 
why the fluctuations 6x are the largest in the latter case 
and turn out to be system dependent, even if the non- 



interacting geometries all exhibit chaotic dynamics. 
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FIG. 3. Unoccupied SCHF levels e)i, ' ,n = N +1,N + 2, . . . 
as a function of interaction strength for a chaotic quantum 
dot (bottom right inset in Fig. El) with A*' =7, 11, 16, and 
24 electrons. Note the sudden jumps in the level spac- 
ings E2 between the lowest two adjacent levels which are 
caused by abrupt changes in the mean-field potential. For 
clarity, the levels in each panel are shifted by an amount 
a{N) + h{N)R/aB. The energy units are h^/2m*R^. 



(iii) The shape dependence of the capacitive energy 
fluctuations found within our model, which does not in- 
clude effects, e.g. from the dot surroundings, tempera- 
ture, and spin, may be considered as one possible mech- 
anism which can lead to the different widths in the peak 
spacing distributions of the various experiments [6-9] . 

An experimental observation of the interaction depen- 
dent crossover of Sx would be particularly interesting. 
For GaAs samples with a^ « 10 nm, the cross-over may 
appear for a system size close to 0.1/im. Recent experi- 
ments |23] on small quantum dots without spatial sym- 
metry may be a good testing ground. 

To conclude, we considered interaction effects in clean 
chaotic quantum dots. We propose as a mechanism for 
enhanced capacitive energy fluctuations the formation of 
ordered charge states near the quantum dot edges. In the 
regime where the ordered states are developed the fluctu- 
ations are geometry-dependent, even for chaotic systems. 

We thank R. Berkovits, H. Castella, A. Cohen, Y. 
Gefen, S. Kettemann, A.D. Mirlin, G. Montambaux, and 
J.-L. Pichard for helpful discussions. KHA thanks M. 
Takahashi for useful comments on numerical calculations. 
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FIG. 4. Charge reordering with increasing interaction for 
a chaotic dot (bottom inset in Fig. 0); R/a% = 5 (15) in 
left (right) panels. Panels (a) and (b) show the probabil- 
ity density of the first unoccupied SCHF single-particle state 
n = 8, N = 7 (lowest curve in upper left panel. Fig. m). In (c) 
to (f) the total charge densities for A'' = 7, 11 are displayed. 

(ii) The accumulation of charge near the edge alone is 
not sufficient to produce enhanced fluctuations Sx, as vis- 
ible in Fig. ^(c, e) for weakly interacting dots {R/ag=5). 
This weak interacting case, where 5x is still similar to the 
RMT value, can be understood within the RPA approach 
of Ref . ||l^ , where the corresponding contribution to Sx 
was shown to be of order of A or smaller, independent of 
the interaction strength. The formation of ordered states, 
which no longer follow RMT predicitions, is necessary to 
generate enhanced fluctuations. 
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